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Abstract. 

We demonstrate an alternative correlation between the El Nino Southern Oscillation 
(ENSO) and global temperature variation to that shown by McLean et al. [2009]. We 
show 50% of the variation in RATPAC-A tropospheric temperature (and 54% of Had- 
CRUT3) is explained by a novel cumulative Southern Oscillation Index (cSOI) term in 
a simple linear regression. We review evidence from physical and statistical research in 
support of the hypothesis that accumulation of the effects of ENSO can produce nat- 
ural multi-decadal warming trends. Although it is not possible to reliably determine the 
relative contribution of anthropogenic forcing and SOI accumulation from multiple re- 
gression models due to collinearity, these analyses suggest that an accumulation ratio 
cSOI/SOI of 4.8 ± 1.5% and up to 9 ± 2% is sufficient for ENSO to play a large part 
in the global mean temperature trend since 1960. 



1. Introduction 

McLean et al. [2009] claim a high correlation between the 
Southern Oscillation Index (SOI) and tropospheric temper- 
ature on lagged, detrended data. These correlations sug- 
gest the involvement of the El Nino Southern Oscillation 
(ENSO) in global temperature variations, but not warming 
trends due to the various detrending operations performed 
in the course of analysis. Foster et al. [2009] question their 
approach, also citing research to date indicating that SOI 
explains between 15 and 30% of interannual variability in 
surface and/or lower tropospheric temperature, but little 
of the global mean warming trend, generally attributed to 
an increase in anthropogenic greenhouse gas concentrations 
[IPCC, 2007]. We bring attention to another suggestion by 
Bob Tisdale (http://bobtisdale.blogspot.com/) for explain- 
ing warming trends via a cumulative function of the SOI 
(cSOI), representing a residual effect of ENSO events on 
global temperature, with potential to account for a large 
part of the increase in global temperature since 1960. 

The increase in the frequency and intensity of ENSO since 
the mid-1970s compared with the previous 100 years may 
have increased global surface temperature during the last 
few decades, although the magnitude of the effect is thought 
to be small [IPCC, 2001]. Recent results show a radiation 
imbalance related to the Pacific Decadal Oscillation (PDO) 
of +0.15W/m 2 over the period [Douglass and Knox, 2009]. 
Various lines of research provide physical explanations for 
globally persistent ENSO-effects [Jacobs et al., 1994; White 
et al, 2003]. In particular, the El Nino events produce larger 
sea surface anomalies more than La Nina [Monahan and 
Dai, 2004] and increase solar absorption [Pavlakis et al., 
2008]. El Nino effects may be due intensified by positive 
feedbacks such as relaxing of the trade winds further increas- 
ing heating [Sun and Yu, 2009], or strong positive feedback 
from low level clouds over the Pacific [Clement et al., 2009]. 
Sun and Yu [2009] propose a dynamic whereby the asym- 
metric warming leads to a slow increase in the Pacific mean 
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temperature during periods of persistent El Nino, ultimately 
leading to a reversion back to more neutral La Nina condi- 
tions. Thus, it has been speculated that global warming on 
decadal, interdecadal and even centennial scales may be due 
to internal modes of variability [White et al., 2003]. The 
potential for a natural source of warming during the 1960- 
1998 period is supported by records of decreasing cloud cover 
during the El Nino dominated positive phase of the Pacific 
Decadal Oscillation [Wielicki et al, 2002; Palle et al, 2005]. 

Purely statistical analysis suggests a (nearly) unit-root 
(or random walk) property in the global temperature series 
[Stockwell, 2006; Koutsoyiannis and Cohn, 2008]. It is well 
known that a random walk can be produced simply from the 
cumulative sum of an independent random variable, and al- 
though the source has not been identified, physically realis- 
tic versions of random walks, such as mean-reverting series, 
long-term persistent (LTP), or self-similar series, are sug- 
gested by global temperature trends [Cohn and Lins, 2005]. 
The apparent unit-root property of global temperature may 
be caused ocean integration of ENSO effects. 

Some climate models show multi-decadal natural varia- 
tion in temperature, although attribution of these trends to 
ENSO may be confounded, as some GCMs show an ENSO- 
like response to global warming (although others show an 
enhanced Northern Annular Mode) [IPCC, 2007]. Models 
with insufficient variability at decadal time scales are un- 
likely to show ENSO accumulation. For example, a mul- 
tilayer ocean model would be needed to represent a persis- 
tent re-emergence mechanism as proposed by Alexander and 
Deser [1995], whereby anomalies carry over from the previ- 
ous years winter water through entrainment into the mixed 
layer as it deepens the following winter [Timlin et al., 2002]. 
Clouds are still acknowledged as the major source of uncer- 
tainty, and inadequacies in the tropical energy budget such 
as the dynamics and variability of top-of-atmosphere radi- 
ation [Wielicki et al, 2002] and lack of realism in the han- 
dling of the tropical precipitation relations and transfer of 
heat [Jungclaus et al, 2006] might also impede recognition. 
Thus while models are broadly consistent in attributing re- 
cent temperature trends to anthropogenic forcing, they do 
not rule out the potential for internal variability to account 
for a large part of the increase in temperature from 1960 to 
2008 [Wigley and Paper, 1995]. 

2. Analysis 
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cSOI(t) = SO h 
t=i 

We fit a series of multiple regression models to 12-month 
filtered (or smoothed) tropospheric global temperature se- 
ries as per McLean et al. [2009] (RATPAC-A available at 
http://wwwl.ncdc.noaa.gov from NOAA [2009]), and to 
an unfiltered global land and ocean surface dataset (Had- 
CRUT3 - available at http://www.cru.uea.ac.uk by Jones 
et al. [1999]). The independent terms in the models were 
composed of the Southern Oscillation Index (SOI - available 
at http://www.bom.gov.au/climate/current/soihtml.shtml), 
the sunspot number as a proxy for solar irradiance (SSN 
- available at http://solarscience.msfc.nasa.gov) and from 
the GISS [2009] website, the stratospheric aerosols (V - 
available at http://data.giss.nasa.gov) and anthropogenic 
forcing from the smoothed sum of all variables other 
than volcanism and solar forcing (AF - available at 
http://data.giss.nasa.gov). The SOI was not lagged. 

Figure 1 shows various fitted models for RATPAC-A (up- 
per, filtered as in McLean et al. [2009]) and HadCRUT3 
(lower). The fit of cSOI (blue line) to temperature is ap- 
parent, and further enhanced with SOI, volcanics (V) and 
sunspot number (SSN) (red line). Inclusion of the anthro- 
pogenic forcing term (AF) increases the overall fit slightly 
(green line). Also worthy of note, the maximum and min- 
imum of global temperatures (HadCRUT3) are coincident 
with maximum and minimum of cSOI over the same period. 

Table 1 lists the correlation results for the multiple regres- 
sion models of global temperature. The R coefficients of 
determination (proportion of variation explained) for the in- 
sample (R2) and the split-segment cross-validation average 
(R2c, split on 1984) are listed for both RATPAC-A (RAT) 
(SOI and RATPAC 12-month filtered but not lagged), and 
HadCRUT3 (HAD) (not filtered or lagged) global surface 
temperature data sets. 

Referring to the results for RATPAC in Table 1, model 
GTTA ~ SOI, though significant, explains a small propor- 
tion of the variation (7%) from 1960 to 2008. The addi- 
tion of the cSOI term increases the variance explained to 
50% (GTTA ~ cSOI). The anthropogenic forcing model, 
GTTA ~ AF, explains 69% of variation, increasing to 71% 
when SOI terms are included. The full model explains 78% 
of the variation. 

The variation explained in the HadCRUT3 data (Table 1) 
by cSOI is a larger 54%, and the inclusion of SOI terms im- 
prove the anthropogenic forcing model from 66% to 72%. 
Finally, the full model explains 77% of the variation. The 
cross-validation R 2 results are lower but consistent with the 
in-sample validation results, demonstrating the models are 
robust to out-of-sample split-segment validation. 

While one would like to know the relative contribution 
of anthropogenic forcing and SOI accumulation, reliable de- 
termination with multiple regression models is not possible 
due to the collinearity of the variables. We can, however, 
calculate some values suggested by these analyses and with 
greater confidence estimate an upper bound on ENSO ac- 
cumulated in the absence of anthropogenic forcing. The 
coefficients for two equations (all very significant) were: 

Ti =0.29 ±0.01 

-0.00055 ± O.OOOOlcSOI - 0.0061 ± 0.0006SOJ 

-1.8 ± 0.2U - 0.0011 ± 0.0001SSW 
T 2 =-0.15 ±0.03 

-0.00025 ± 0.00002c5O7 - 0.0051 ± 0.0005SO/ 

-1.6 ± 0.2U - 0.0007 ± 0.000 ISSN + 0.19 ± 0.01AF 



The accumulating proportion of SOI is estimated from 
the ratio of the coefficients of cSOI and SOI, and from the 
maximum ranges of temperature and cSOI. The ratio of co- 
efficients cSOI/SOI for HadCRUT3 (Ti) indicates that up 
to 9±2% (95%CL from coefficient uncertainty) of the SOI is 
accumulated. The equivalent ratio of coefficients cSOI/SOI 
for HadCRUT3 (Tb) including anthropogenic forcing (AF) 
suggests a value of 4.9 ± 1.5% may represent the fraction 
of SOI accumulated, however this figure needs to be viewed 
with caution. 

The cSOI ranged from 948 in 1976 to -172 in 1998 and 
global temperature ranged 0.8817 degrees Celsius between 
the coincident maximum and minimum over the same pe- 
riod. Therefore the contribution of cSOI to the increase is 
potentially 0.8817/1120 or 0.00079 degrees of global tem- 
perature per SOI unit. When divided by the SOI coefficient 
(Ti) this gives a higher estimate of 13 ± 3%, as it incor- 
porates some 'overshoot' at the extreme highs and lows of 
temperature. 

3. Discussion 

These analyses suggest that a contribution from ENSO- 
effects to global temperatures, when expressed as the cumu- 
lative sum of the SOI, can potentially account for 50% of the 
variation in global mean temperature in the last 50 years - a 
'large part' of warming, as claimed by McLean et al. [2009]. 

We determine the relative contribution of ENSO and an- 
thropogenic forcing, with the qualification that linear regres- 
sion studies are not reliable when variables are collinear, and 
other data are needed to determine whether both, either, or 
neither contribute in reality. However, our estimated sensi- 
tivity of global mean temperature to anthropogenic forcing 
when cSOI is included is 0.19 ± 0.02%K/W/m 2 (AF coeffi- 
cient equation 2), which is smaller, though close to the no- 
feedback sensitivity value for a number of forcings (including 
CO2) of 0.30K/W/m 2 [Douglass and Christy, 2008]. As an 
estimated sensitivity of around 0.8K/W/m 2 is required for 
climate models to reproduce the observed warming trend 
since 1960, it follows that if the accumulated ENSO effect 
were to make up the difference in global mean trend since 
1960, ENSO must be responsible for greater than 50% of 
the observed warming. Therefore, providing interactions be- 
tween ENSO and greenhouse gas forcing are not significant, 
these results are clearly inconsistent with the claim that 
" [M]ost of the observed increase in global average tempera- 
tures since the mid-twentieth century is most likely due to 
the observed increase in anthropogenic greenhouse gas con- 
centrations" [IPCC, 2007]. The global temperature change 
from ENSO events that needs to be accumulated to achieve 
the observed trends is estimated from the ratio of the coef- 
ficients cSOI/SOI at 4.9 ± 1.5% and may be up to 9%. 

These analyses show it would be premature to dismiss a 
large role for ENSO in the global mean temperature trend 
since 1960. The cumulative SOI together with physical re- 
search and statistical modelling, provide support for the con- 
tention that multi-decadal natural variability from ENSO 
events may play a larger part in recent warming than is gen- 
erally appreciated, although further studies of longer tem- 
poral duration and spatial extent are needed to better un- 
derstand and quantify the mechanisms involved. 
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Table 1. In-samplc (R2) and split-sample cross-validation 
(R2c) correlation coefficients of regressions of global mean 
temperature series RATPAC-A (RAT 12-month filtered) and 
HaDCRUT3 (HAD). 
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Figure 1. The fit of models to temperature for 
RATPAC-A filtered (above) and HadCRUT3 raw (be- 
low), including the global temperature (gray line), cumu- 
lative SOI only (blue line), the natural model (red line), 
and full model including anthropogenic forcing (green 
line) . 



